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'fi\" We present an improved calculation of all B — > light pseudoscalar formfactors 

f-H ' from light-cone sum rules, including one-loop radiative corrections to twist-2 and 

twist-3 contributions, and leading order twist-4 corrections. The total theoretical 
uncertainty of our results at zero momentum transfer is 10 to 13%. The dependence 
1 J - of the formfactors on the momentum transfer g^ is parametrized in a simple way 

|^N( ' that is consistent with their analytical properties and is valid for all physical g^. 

'H The uncertainty of the extrapolation in q^ on the semileptonic decay rate T{B ^ 

ireu) is estimated to be 5%. 



1. Introduction and Definitions 

In a recent paper^ we have reported a new calculation oi B ^ ir, K,rj decay 
formfactors from QCD sum rules on the light-cone (LCSRs). The paper 
improves upon our previous publications^''^ by: 

• including radiative corrections to twist-3 contributions to one-loop 
accuracy, for all formfactors; 

• a precisely defined method for fixing sum rule specific parameters; 

• using updated values for input parameters; 

• a careful analysis of the uncertainties of the formfactors at zero 
momentum transfer; 

• a new parametrization of the dependence of the formfactors on 
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monientuni transfer, which is consistent with the constraints from 
analyticity and heavy-quark expansion; 
• a detailed breakdown of the dependence of formfactors on non- 
perturbative hadronic parameters describing the tt, K^ t] mesons, 
which facihtates the incorporation of future updates of their nu- 
merical values and also allows a consistent treatment of their effect 
on nonleptonic decays. 

The key idea of LCSRs is to consider a correlation function of the weak cur- 
rent and a current with the quantum-numbers of the B meson, sandwiched 
between the vacuum and, in the present context, the pseudoscalar meson P, 
i.e. TT, K and rj. For large (negative) virtualities of these currents, the corre- 
lation function is, in coordinate-space, dominated by distances close to the 
light-cone and can be discussed in the framework of light-cone expansion. 
In contrast to the short-distance expansion employed by conventional QCD 
sum rules a la SVZ^, where nonperturbative effects are encoded in vacuum 
expectation values of local operators with vacuum quantum numbers, the 
condensates, LCSRs rely on the factorisation of the underlying correlation 
function into genuinely nonperturbative and universal hadron distribution 
amplitudes (DAs) (j) that are convoluted with process-dependent amplitudes 
T, which arc the analogues to the Wilson-coefficients in the short-distance 
expansion and can be calculated in perturbation theory. Schematically, one 
has 

correlation function - ^ T^") (g) 0("). (1) 

n 

The sum runs over contributions with increasing twist, labelled by n, which 
are suppressed by increasing powers of, roughly speaking, the virtualities 
of the involved currents. The same correlation function can, on the other 
hand, be written as a dispersion-relation, in the virtuality of the current 
coupling to the B meson. Equating dispersion-representation and the light- 
cone expansion, and separating the B meson contribution from that of 
higher one- and multi-particle states, one obtains a relation (QCD sum 
rule) for the formfactor describing the B ^ P transition. 

The particular strength of LCSRs lies in the fact that they allow the 
inclusion not only of hard-gluon exchange contributions, which have been 
identified, in the seminal papers that opened the study of hard exclusive 
processes in the framework of perturbative QCD (pQCD)'"', as being dom- 
inant in light-meson form factors, but that they also capture the so-called 
Feynman-mechanism, where the quark created at the weak vertex carries 
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nearly all momentum of the meson in the final state, while all other quarks 
are soft. This mechanism is suppressed by two powers of momentum- 
transfer in processes with light mesons, but there is no suppression in 
heavy-to- light transitions®, and hence any reasonable application of pQCD 
to B meson decays should include this mechanism. It is precisely LCSRs 
that accomplish this task and have been applied to a variety of problems 
in heavy-meson physics. ^'^'^^^ A more detailed discussion of the rationale 
of LCSRs and of the more technical aspects of the method can be found 
e.g. in Ref.*. 

The formfactors in question can be defined as {q = pB — p) 

{P{p)\qj,b\B{pB)) = f^iq') |(PB +p), - ^^^^^^ Q, 

^' .foil')l,, (2) 



q' 



{P{p)\q<J^uq''b\B{pB)) =1{{PB +P)^q^ - q^{ml - m%)} '^'^ ' . (3) 

TUB + tnp 

The starting point for the calculation of e.g. /" is the correlation function 

»y dV^^(7r(p)|r[q-7^6](y)K&z75<z](0)|0) = n+2p^ + . . . , (4) 

where the dots stand for other Lorentz structures. For a certain config- 
uration of virtualities, namely m^ — P% > 0{AQCDnib) and ml — q^ > 
0(Aqcd"t-6), the integral is dominated by light- like distances and can be 
expanded around the light-cone: 

n+{q\pl)=J2f rf^i0(")(w;A*F)T(")(w;g^p|;MF). (5) 

As in Eq. (1), n labels the twist of operators and /ip denotes the factori- 
sation scale. The restriction on g^, ml — q^ > ©(AqcD^f,), imphes that 
/^ is not accessible at all momentum-transfers; to be specific, we restrict 
ourselves to < g^ < 14GeV^. As 11+ is independent of y^p, the above 
formula implies that the scale-dependence of T*^"^ must be canceled by that 
of thcDAs 0("). 

In Eq. (5) it is assumed that 11+ can be described by coUinear factori- 
sation, i.e. that the only relevant degrees of freedom are the longitudinal 
momentum fractions u carried by the partons in the tt, and that transverse 
momenta can be integrated over. Hard infrared (coUinear) divergences oc- 
curring in r'"-* should be absorbable into the DAs. CoUinear factorisation is 
trivial at tree-level, where the h quark mass acts effectively as regulator, but 
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can, in principle, be violated by radiative corrections, by so-called "soft" 
divergent terms, which yield divergences upon integration over u. Actually, 
however, it turns out that for all formfactors calculated in Ref.^ the T are 
nonsingular at the endpoints u = 0, 1, so there are no soft divergences, in- 
dependent of the end-point behavior of the distribution amplitudes. In Ref.^ 
Eq. (5) has been demonstrated to be valid to 0{as) accuracy for twist-2 
and twist-3 contributions for all correlation functions n_|_^o,T from which to 
determine the formfactors /+.o,t- 

As for the distribution amplitudes (DAs), they have been discussed 
intensively in the literature^. For pseudoscalar mesons, there is only one 
DA of leading- twist, i.e. twist-2, which is defined by the following light-cone 
matrix clement (a;^ = 0): 

(0|u(x)7^75rf(-2;)k(p))-*./.P^ / d7ie^^P^0^(u), (6) 

Jo 

where (^ = 2u — 1 and wc have suppressed the Wilson-line [x, —x] needed 
to ensure gaugc-invariance. The higher-twist DAs are of type 

(0\u{x)rd{~x)\TT{p)) or (0|M(x)rG°^(wa;)A"/2d(-.T)|7r(p)), 

where w is a number between and 1 and T a combination of Dirac matrices. 
The sum rule calculations performed in Refs.^'^^'^ include all contributions 
from DAs up to twist-4. The DAs are parametrized by their partial wave 
expansion in conformal spin, which to NLO provides a controlled and eco- 
nomic expansion in terms of only a few hadronic parameters^. 

The LCSR for /^ is derived in the following way: the correlation func- 
tion n_|_ , calculated for unphysical p^ , can be written as dispersion relation 
over its physical cut. Singling out the contribution of the B meson, one has 

n+ = /I(g^) 7^ o + higher poles and cuts, (7) 

^b-Pb 

where /b is the leptonic decay constant of the B meson, fB''n% = 
mb{B\bijzd\0). In the framework of LCSRs one docs not use (7) as it 
stands, but performs a Borel transformation, l/(t — p^) -^ B l/{t — p^) = 
1/M^exp(— i/M^), with the Borel parameter M^; this transformation en- 
hances the ground-state B meson contribution to the dispersion represen- 
tation of n+ and suppresses contributions of higher twist to the light-cone 
expansion of 11+. The next step is to invoke quark-hadron duality to ap- 
proximate the contributions of hadrons other than the ground-state B me- 
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son by the imaginary part of the light-cone expansion of 11+ , so that 

1 r^ 

- dtlmnLCE(t) cxp(-t/M2). (8) 



Af2 

1 1 



Subtracting the 2nd term on the right-hand side from both sides, one ob- 
tains 

Eq. (9) is the LCSR for /^. sq is the so-cahed continuum threshold, which 
separates the ground-state from the continuum contribution. At tree-level, 
the continuum-subtraction in (9) introduces a lower limit of integration in 
u, the momentum fraction of the quark in the tt: u > (ml — q'^)/{so — q^), 
in (5), which behaves as 1 — Aqcd /'^it, for large mi, and thus corresponds 
to the dynamical configuration of the Feynman-mechanism, as it cuts off 
low momenta of the u quark created at the weak vertex. At 0{as), there 
are also contributions with no cut in the integration over u, which corre- 
spond to hard-gluon exchange contributions. The task now is to find sets 
of parameters M^ (the Borel parameter) and so (the continuum threshold) 
such that the resulting formfactor does not depend too much on the precise 
values of these parameters. 



2. Results 

For a detailed discussion of the procedure used to determine the hadronic 
and sum rule specific input parameters we refer to Ref.^. One main fea- 
ture is that /b, the decay constant of the B meson entering Eq. (9) is 
calculated from a sum rule itself^", which reduces the dependence of the 
resulting formfactors on the input parameters, in particular vrih, which is 
the one-loop pole mass and taken to be (4.80 ± 0.05) GeV. This procedure 
does not, however, reduce the formfactors' dependence on the parameters 
describing the twist-2 DAs, which turns out to be rather crucial. Despite 
much effort spent on both their calculation from first principles and their 
extraction from experimental data, these so-called Gegenbauer moments, 
fli (only for K), 02 and 04 (for all P) are not known very precisely. Figure 1 
shows the dependence of /^(O) on 02 and 04; the dots represent different 
determinations of these parameters and illustrate the resulting spread in 
values of the formfactor. The situation is even more disadvantageous for 
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Figure 1. Dependence of /7!(0) on a2 and 04, for central values of input parameters. 
The lines are lines of constant /T(0). The dot labeled BZ denotes our preferred values 
of 12, 4i BMS the values from the nonlocal condensate model^^ and BF from sum rule 
calculations^ . 





Figure 2. (a) Dependence of /i!^(0) on the Gegenbauer moment a\. (b) f^{q^) as 
function of q^ for different values of ay. solid line: a^ = 0.17, short dashes: a^ = 0, 
long dashes: a^ = —0.18. 



the if, whose fornifactors depend on the SU(3) breaking parameter af^, 
whose size and even sign are under discussion -'^^: at present, values as dif- 
ferent as —0.18 and +0.17 (at /i = 1 GeV) are being quoted. Figure 2 shows 
the dependencies of (a) f^{0) and (b) f^{q^) on this parameter; evidently 
it is very important to determine its value more precisely. 

Summarizing the detailed analysis of the uncertainties induced by both 
external input and LCSR parameters, the final results for the formfactors 
at zero momentum transfer obtained in Ref.^ are: 



/^(O) = 0.258 ±0.031, 
/f(0) = 0.331 ±0.041 
fliQ) = 0.275 ±0.036, 



0.25(5a,, 



/|(0) = 
/?(0) = 



0.253 ±0.028, 
0.358 ±0.037 
0.285 ±0.029. 



0.31<5a,, 



5ai is defined as af^(l GcV) — 0.17, i.e. the deviation of af from the cen- 
tral value used in Ref.^. For f^-"^ the total theoretical uncertainty ranges 
between 10% to 13%, for /^ it is 12%, plus the uncertainty in ai, which 
hopefully will be clarified through an independent calculation in the not 
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Figure 3. /+ (solid lines), /o (short dashes) and /y (long dashes) as functions of (/^ for 
TT, K and r). The renormalisation scale of fx is chosen to be m;,. 

too far future. The intrinsic, irreducible uncertainty of the sum rule cal- 
culation is related to the dependence of the result on the sum rule specific 
parameters M^ and Sq and estimated to be ^ 7%. 

Turning to the (/^-dependence of formfactors, it has to be recalled that 
LCSRs are only valid if the energy Ep of the final state meson, measured 



in the rest frame of the decaying B, is large, i.e. if q^ 



— 2mBEp is 



not too large; specifically, we choose Ep > 1.3 GeV, i.e. q^ < 14GcV^. The 
resulting formfactors are plotted in Fig. 3, using central values for the input 
parameters. In order to allow a simple implementation of these results in 
actual applications, and also in order to provide predictions for the full 
physical regime < q^ < {itib — mp)'^ « 25 GeV , it is necessary to find 
parametrizations of /(g^) that 

• reproduce the data below 14 GeV with good accuracy; 

• provide an extrapolation to q^ > 14 GeV that is consistent with 
the expected analytical properties of the formfactors and repro- 
duces the lowest-lying resonance (pole) with J^ = 1^ for /+ and 



^For /(), the lowest pole with quantum numbers 0+ lies above the two-particle threshold 
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Table 1. Fit parameters for f{q'^). mi is 
the vector meson mass in the correspond- 
ing channel: m^'^ = nig* = 5.32 GeV and 



= mg, = 
: 4. 8 GeV. 



5.41 GeV. The scale of fr is 





n 


r2 


(mi)2 


<t 


fl 


0.744 


-0.486 


{mlY 


40.73 


f^ 





0.258 


- 


33.81 


f^ 


1.387 


-1.134 


(mj)2 


32.22 


/f 


0.162 


0.173 


(mf)2 


- 


/n^ 





0.330 


- 


37.46 





0.161 


0.198 


(mf)2 


- 


n 


0.122 


0.155 


K)2 


- 


fV 

Jo 





0.273 


- 


31.03 


n 


0.111 


0.175 


K)2 


- 



As shown in Ref.^, the foUowing paranietrizations are appropriate: 



for fl 



,t' 



J \^ ) 1 2/2~'~i 2/2 

1 — q^'/rrii 1 — q'^ /m^ 



(10) 



'fit 



where ml is the mass of B*{1 ), ml ~ 5.32 GeV; the fit parameters 
are ri, r2 and mfit; 

• for /^'i?: 

r( 2^ ^ ^1 ^2 .... 

1 — q-'/mf (1 — q^ /m,\y 

where Wi is the mass of the 1^ meson in the corresponding channel, 
i.e. 5.32 GeV for 77 and 5.41 GcV for K\ the fit parameters are ri 
and r2] 

• for fo- 



r2 



1 - q'/mi^ 
the fit parameters are r2 and mfit. 



(12) 



The central results for the fit parameters are collected in Tab. 1. The 
quality of all fits is very good and the maximum deviation between LCSR 
and fitted result is 2% or better. The impact of the extrapolation of the 
fit formulas to q^ > 14 GeV^ is of phenomenological relevance mainly for 
B -^ TTev, relevant for the determination of \Vub\ from experiment. We have 



starting at (mg + mp)'^ and hence is not expected to feature prominently. 
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estimated the effect of the extrapolation on the decay rate by implementing 
different parametrisations for /^, which all fit the LCSR result very well for 
q^ < 14GeV , but differ for larger g^, the main distinguishing feature being 
the positions of the poles. We find that for reasonable parametrisations, 
that is such that do not exhibit too strong a singularity at q^ = m^, the 
total rates differ by not more than 5%, the difference becoming smaller if an 
cut-off on the maximum invariant mass of the lepton pair is implemented, 
which implies that the extrapolation is well under control. 



3. Summary & Conclusions 

LCSRs provide accurate results for weak decay formfactors of the B meson 
into light mesons, in particular tt, K and rj. The results depend on sum 
rule specific input parameters which generate an irreducible "systematic" 
uncertainty of the approach estimated to be '--^ 7%. Additional uncertainties 
are induced by imprecisely known hadronic input parameters, in particular 
the Gcgenbaucr moments ai.2,4 describing the leading-twist light-meson 
distribution amplitudes. An improved determination of these parameters 
would be very welcome. The present total uncertainty of the formfactors at 
zero momentum transfer varies between 10 and 13%, but becomes smaller 
at larger q^. LCSR calculations require the energy of the final state meson 
to be large in the rest-frame of the decaying B and hence are valid only for 
not too large momentum transfer q^; the maximum eligible q^ is chosen to 
be 14 GeV^ . The g^-dependence of the formfactors can be cast into simple 
parametrizations in terms of two or three parameters, which also capture 
the main features of the analytical structure and are expected to be valid 
in the full kinematical regime < g^ < [niB — mp)^ . The total uncertainty 
introduced by the extrapolation of the formfactors to q^ larger than the 
sum rule cut-off 14 GeV^ is estimated to be ~ 5% for the semileptonic rate 
T{B -^ TTev). 

Ref.^ also contains a detailed breakdown of the dependence of the form- 
factors on the Gegenbauer moments, which not only allows one to recalcu- 
late the formfactors once these parameters are determined more precisely, 
but also makes it possible to consistently assess their impact on nonlcptonic 
decay amplitudes (e.g. B -^ tttt) treated in QCD factorisation. 

The LCSR approach is complementary to standard lattice calculations, 
in the sense that it works best for large energies of the final state meson 
(i.e. small q^), whereas lattice calculations work best for small energies - a 
situation that may change in the future with the implementation of moving 
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NRQCD^'^. Previously, the LCSR results for /^ q at small and moderate q^ 
were found to nicely match^^ the lattice results obtained for large q^ . The 
situation will have to be reassessed in view of our new results and it will 
be very interesting to see if and how it will develop with further progress 
in both lattice and LCSR calculations. 
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